Quantum breaking time near classical equilibrium points 
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By using numerical and semiclassical methods, we evaluate the quantum breaking, or Ehrenfest 
time for a wave packet localized around classical equilibrium points of autonomous one-dimensional 
systems with polynomial potentials. We find that the Ehrenfest time diverges logarithmically with 
the inverse of the Planck constant whenever the equilibrium point is exponentially unstable. For 
stable equilibrium points, we have a power law divergence with exponent determined by the degree 
of the potential near the equilibrium point. 
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The question of estimating how long classical and 
quantum evolutions stay close is one of the main prob- 
lems of semiclassical analysis. The evolution of a quan- 
tum observable can follow that of the corresponding clas- 
sical one up to a finite time, the so called quantum break- 
ing, or Ehrenfest, time. As initially conjectured in fj], |^| 
and rigorously proved in ||, ||, |5| , whenever the classical 
flow is chaotic, the Ehrenfest time diverges logarithmi- 
cally in h. This result is easily understood. Starting from 
an initial value A(ft) ~ ft//, where / is a characteristic ac- 
tion of the system, the difference between a classical flow, 
with Lyapunov exponent A > 0, and the corresponding 
quasi-periodic quantum flow increases as A(fi) exp(At). 
The two flows depart at t ~ A -1 \og(I/h). The situation 
is different for a regular classical flow. In this case, start- 
ing from the work ||, it was suggested in that the 
Ehrenfest time grows algebraically as ftr & . The determi- 
nation of the value of 5 and its possible universal nature 
is still an open problem. See J^] and references therein 
for recent results. 

The ft-scaling of the Ehrenfest time is usually investi- 
gated for classical flows which are completely chaotic or 
regular. However, it is interesting to study the quantum- 
classical correspondence in systems having isolated un- 
stable orbits embedded in a completely regular phase- 
space. The simplest example is given by the ubiquitous 
double- well system defined by the Hamiltonian H(p, q) — 
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2 + • For this system, there is only one unstable 



periodic orbit, namely that associated to the equilibrium 
point (po:9o) = (0,0), with positive Lyapunov exponent 
A = 1. Is it possible to have a logarithmic Ehrenfest time 
in proximity of an isolated exponentially unstable point 
like (po,?o)? 

The usual way of studying the Ehrenfest time consists 
in comparing the evolution of classical observables with 
the quantum expectation value of the corresponding op- 
erators, either in the coherent state representation || , or 
in the framework of Weyl quantization . In the present 
case, we follow a simpler approach based on the analy- 
sis of the quantum spectrum. We know that on going 



towards the classical equilibrium point (po,9o) the pe- 
riod of motion diverges, so that the evolution of a phase- 
space distribution function localized around this point 
must show a continuous frequency distribution around 
v = 0. On the other hand, in the quantum case, due to 
the discrete nature of the spectrum, the frequency dis- 
tribution is characterized by a gap between zero and a 
minimal frequency. We call this minimal frequency the 
Ehrenfest frequency, ve- In fact, its inverse, v^ 1 , is an 
upper bound to the time at which the quantum-classical 
correspondence of the evolution of any observable breaks 
down. We estimate the Ehrenfest time as . 

By using numerical and semiclassical methods, we 
study the behavior of ve{K) around classical equilibrium 
points, both stable and unstable, for several autonomous 
one-dimensional systems. We find that v~j^(K) diverges 
logarithmically for h — > whenever the equilibrium point 
is exponentially unstable. In all the other cases, the 
Ehrenfest time follows a power law with exponent re- 
lated to the degree of the potential near the equilibrium 
point. 

In the following, we consider systems described by the 
Hamiltonians 
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H{ V ,q) = ?-+A q —+B?-, 

2m 2a 2p 



(1) 



with A < 0, B > and (3 > a > 1. By properly rescaling 
position, momentum and time, we can always reduce to 
the case B = 1, m = 1 and either A = or A = — 1 |fl6| . 
For A — 0, we have single- well systems with a classical 
stable equilibrium point (po, qo) = (0,0) at energy e = 0. 
A more interesting situation occurs for A = — 1. In this 
case the systems are double- well oscillators and the clas- 
sical equilibrium point (po, qo) = (0,0) at energy e = is 
unstable. In the particular case a = 1, the equilibrium 
point is exponentially unstable. In both cases, A = or 
A = — 1, the periodic orbits near the equilibrium point 
at e = have a period which diverges for e — » 0. 

On the quantum mechanical side, in order to repre- 
sent a state localized near the classical equilibrium point 
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FIG. 1: Fourier transform of the survival probability V(v) for 
different values of h in the case A = 0, j3 = 2. The dashed 
line is the h — * limit distribution given by Eq. (jl4|) . 



(poj <7o) we consider the following initial wavefunction 



(«l^(0)) 
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The associated Wigner function, 
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Wfy(P)9) = -rexp 
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(3) 



can be interpreted as a phase-space probability distri- 
bution centered around the point (pojffo) and has the 
property 



lim W^,(p,q) = S(p-p )S(q - q ). 



(4) 



In these expressions fi. is the adimensional rescalcd Planck 
constant, which vanishes when, for instance, the mass m 
of the system is taken larger and larger. 

Instead of considering the evolution of a specific ob- 
servable, we study the simpler survival probability 



v(t) = \(m\m)f 



(5) 



which contains the same gross dynamical information. In 
the basis of the eigenstates of the Hamiltonian, 



H\cj> n ) = e n \(j) n ), n = 0,l,2,..., 
the survival probability V(t) can be written as 



= E E i c « 



I exp (iu nm t) , 



(6) 



(7) 



n— m— 



where c fe = (V>(O)|0k) and v nm = (e„ - £ m )/(2irh). Note 
that Cfe = for k odd, due to the symmetry of the system 
and of the initial wavefunction. 



2 - 



e 4 2 = 0.00 08 148 




2 - 



£40 



-0.0054148 




FIG. 2: Eigenfunctions 040 and 042 corresponding to the min- 
imal frequency ve in the double- well case a = 1, /3 = 2 for 
7i = 10" 2 . The dashed curve is the initial wavefunction (fcj) 
with (po,<7o) = (0,0). 



By using semiclassical and numerical techniques, we 
now show that the Fourier transform of the survival prob- 
ability, 



V(v) = \Cn\ 2 \Crn\ 2 5{v - V nm ), 



(8) 



n— m— 



for sufficiently small values of h is characterized by a gap, 
large with respect to the typical level spacing, between 
v = and a frequency which we call the Ehrenfest fre- 
quency, defined as 



ve 



nun 

ldl 2 l<=ml 2 *0 



(9) 



In the simple case A — 0, by using standard WKB 
approximations, we have 



n+i)/M( /i 



2/3 



with 



*C0 = t/o 



r i 



and 



F( 2^ ) sin a{e.h;(3) 



(10) 



(11) 



(12) 



with a(e,h;l3) = 2^/2(2(3)^ h^e^ . The behavior of 
■p(f ) obtained by using these expressions for e„ and |c £n | 2 
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FIG. 4: Fourier transform of the survival probability V{v) for 
FIG. 3: Superposition coefficients |c n | 2 as a function of £ n /h different values of h in the double-well case a = 1, /3 = 2. 
for h = 10 -3 in the double-well cases a — 1, = 2 (x), a = 2, 
/3 = 4 (□), and a = 3, /? = 6 (O). 



is shown in Fig. [j] in the case (3 = 2. We see that for ft — > 
the frequency distribution 'P(i') approaches a continuous 
limit given by 

P (u) = lim / d £l d£ 2 P (e 1 )p(e 2 )s(u~ £ -^^-) , (13) 

where p(e) = \ c e | 2 4^ and n(s) is obtained by inverting 
e =£„. By using (fl0|) and fll^), we find 

7>oM - 4^0(4^1^1), (14) 

where K is the Bessel function of zero-th order. Fig- 
ure |l| also shows the presence of the gap at v = and its 
shrinking as ft — ► 0. Since the level spacing e n +i — e n in- 
creases by increasing n, the Ehrcnfest frequency (|9|) t urns 
out to be ve = (£2 — £o)/(27rft). According to (|l0|), its 
inverse diverges as 

"e 1 -^. (15) 

We now consider double- well systems, i. e., the case 
A = —1. For these systems, the standard WKB ap- 
proximation fails near the unstable equilibrium point at 
energy e = 0. Only in the particular case a = 1, a reg- 
ularized semiclassical approximation has been developed 
||, |l^, [ll} and the quantization condition for the energy 
levels reads 

, 1 = = cos(^(e, ft)), (16) 
.^1 + exp^f 

where 

«'-»>-ffl-s h|, s-"« r G +i s)-'- <17) 



For the associated eigenfunctions only microlocal expres- 
sions are available []l0| which do not allow for a direct 
determination of the superposition coefficients |c £ | 2 . For 
this reason, in all cases a > 1 we determine numerically 
the eigenvalues and eigenfunctions of the system. With 
standard numerical techniques, this represents an unsur- 
mountable task since the interesting eigenstates, namely 
those close to energy e = 0, have a quantum number n 
which diverges quickly for h — > 0. We bypass the prob- 
lem by using the algorithm which allows to evalu- 
ate selected eigenstates having a very large number of 
nodes. In Fig. || we show, as an example, the couple of 
even eigenfunctions with energy closest to e = in the 
double- well case a = 1, (3 = 2 evaluated for h = 10~ 2 . 
Note that, already for this still relatively large value of 
ft, the corresponding quantum number is n ~ 40. In our 
numerical calculations we go beyond n ~ 10 . 

In Fig. we show the superposition coefficients eval- 
uated for different double- well systems for ft = 10~ 3 . 
We see that |c„| 2 decreases exponentially departing from 
e = 0. For smaller values of ft, the superposition coef- 
ficients |c ra | 2 follow approximately the same exponential 
behavior as a function of \e n \/h and become denser and 
denser. 

The Fourier transform of the survival probability (@) is 
determined by using the eigenvalues and the superposi- 
tion coefficients obtained numerically. In Fig. ^ we show 
V(y) in the case a = 1, j3 = 2 for different values of 
ft. As in the single-well case, at v = we have a gap 
whose width shrinks as ft — > 0. The width of this gap, 
namely the Ehrenfest frequency, is yielded by a couple of 
even consecutive eigenvalues, close to the energy e = of 
the classical equilibrium point. This can be understood 
roughly in the following way. Consider the number of 
states, A4 , in the energy range [e — ft, e + ft] . The frequen- 
cies associated to the eigenvalues in this energy range are 
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spectrum in turn coincide with the energies of the bound 
states of the corresponding confining inverted potential. 

In conclusion, we have shown that the presence of iso- 
lated exponentially unstable orbits is sufficient to break 
the quantum-classical correspondence at a time scale log- 
arithmic in hr 1 . This feature may be relevant in all meso- 
scopic systems which are modeled by one-dimensional 
multi-well Hamiltonians [|l3[ In these systems the 

Ehrcnfcst time behavior is related to experimentally de- 
tectable properties as the classical to quantum crossover 
of the shot noise p5[ . 

We would like to thank Thierry Paul for very stimulat- 
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Cofinanziamento MURST protocollo MM02263577JQ01. 



FIG. 5: Inverse of the Ehrenfest frequency, function 
of h in the double-well cases a = 1, /3 = 2 (x), a — 2, — 4 
(□), and a = 3, j3 — 6 (O). The solid line is the regularized 
WKB prediction based on ( pi||l7| ), while the dashed and dot- 
dashed lines are numerical fits. 



v ~ N~ , so that, in the limit h — > 0, v vanishes if N e di- 
verges. According to Weyl formula, Af e is proportional to 
the classical phase-space volume bounded by the energy 
shells H(p, q) = e±h. This volume can be evaluated ex- 
actly in terms of simple functions in the single-well case 
and in terms of special functions for double- well systems. 
In all cases, we have that Af e diverges when h — > only 
for £ = 0. In the double-well systems, for a = 1, /3 = 2 
the couple of closest eigenvalues has energies of opposite 
sign, as shown in Fig. |^. For a > 1 these eigenvalues are 
both positive if h is sufficiently small. 

The scaling of v^ 1 with h is shown in Fig. || for dif- 
ferent double-well systems. The plotted points are calcu- 
lated using the numerically determined spectrum while 
the solid line represents the inverse of the Ehrenfest fre- 
quency as determined by using the quantization condi- 
tion @|f7|). The Ehrenfest time increases logarithmi- 
cally with hr 1 only in the case a = 1, /3 = 2, i. e., when 
the equilibrium point is exponentially unstable. In all 
the other cases, a numerical fit suggests that 



(18) 



This is the same scaling law which we would obtain, as 
described by Eq. (|i~5|) , in the case of a single- well poten- 



tial V(q) — q / (2a). This fact can be understood in the 
following way. For h — > 0, the discrete eigenvalues of the 
double- well above e = correspond to the energies of the 
continuous spectrum of the barrier —q 2a /(2a) at which 
the transmission coefficient is maximum. According to 
WKB approximation, these resonances of the continuous 
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